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Abstract: Associated Lame potentials V{x) = a{a + l)msn'^ {x , m) + b{b + l)mcn^(x, m)/dn'^(x, m) are 
used to construct complex, PT-invariant, periodic potentials using the anti-isospectral transformation 
X ^ ix + P, where /? is any nonzero real number. These PT-invariant potentials are dehned by V^'^{x) = 
—V{ix + P), and have a different real period from V{x). They are analytically solvable potentials with a 
finite number of band gaps, when a and b are integers. Explicit expressions for the band edges of some of 
these potentials are given. For the special case of the complex potential U^^(x) = — 2msn^(ia; + f3, m), we 
also analytically obtain the dispersion relation. Additional new, solvable, complex, PT-invariant, periodic 
potentials are obtained by applying the techniques of supersymmetric quantum mechanics. 



In the past few years, Bender and others ^ |2j have looked at several complex potentials with PT- 
symmetry and have shown that the energy eigenvalues are real when PT-symmetry is unbroken, whereas 
they come in complex conjugate pairs when PT-symmetry is spontaneously broken. Recently, Mostafazadeh 
[3] has clarified this issue by showing that there is a more general setting of pseudo-hermiticity (of which 
PT-symmetry forms a special case) in which eigenvalues are either real or occur in complex conjugate pairs. 
Scattering problems with complex, PT- invariant potentials have also been investigated UIIHl- However, 
there have been very few papers discussing periodic potentials with PT-symmetry. Only two types of 
PT-invariant, periodic potentials have been considered in detail in the literature, namely, i sin^^"'"^(x) [0] 
and delta function potentials with complex couplings [71IH]. These potentials have been shown to possess 
real band spectra with an infinite number of band gaps. It might be noted here that obtaining these results 
required extensive numerical analysis and no analytic results for band edge energies were possible |9|. The 
purpose of this letter is to construct and study several new classes of analytically solvable, complex, PT- 
invariant, periodic potentials with a finite number of band gaps. Our approach will consist of (i) making 
use of the anti-isospectral transformation x ^ ix + (5 |lUj . and (ii) constructing supersymmetric partner 
potentials from techniques developed in supersymmetric quantum mechanics 

Anti-Isospectral Transformations: We begin with the simple observation that if ijj{x) is a solution 
of the Schrodinger equation for the real potential V{x) with energy E, then ip{ix + fi) \s a. solution of 
the Schrodinger equation for the complex potential —V{ix + (3) with energy — i?, where (3 is an arbitrary 
constant. The new potential —V{ix -\- (3), generated by the anti-isospectral transformation x — > ix -|- /?, is 
clearly PT-symmetric and will be denoted by V^'^{x). Further, if ijj{x) and ^(zx -|- (3) satisfy appropriate 
boundary conditions, they are eigenfunctions of V{x) and V^'^{x) respectively. Since the ordering of 
energy levels for V^'^{x) is the opposite of the ordering of energy levels for V{x)^ this is presumably why 
the transformation x ^ ix + (3 is called "anti-isospectral" . 

In this letter, we focus on periodic potentials by choosing V{x) to be the associated Lame potential 

V{x) = a{a + l)msn2(x, m) + b(b + 1)^ ™^"^'"^^ , (1) 

dn [x, m) 

which has a real period 2K{m). Note that if either a or 6 is zero, this potential is called the Lame potential. 
Recall that when a and b are non-negative integers, the associated Lame potential has many analytically 
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solvable eigenstates and only a finite number of band gaps ^lEJ. Here, sn(j;,m), cn(x,m), dn(3;,m) 
are Jacobi elliptic functions with elliptic modulus parameter m (0 < m < 1). They are doubly periodic 
functions with periods [4K{m),i2K'{m)], [4K{m),2K{m) + i2K'{m)], [2K{m),i4:K' (m)] respectively 
[imn], where K{m) = de[l — m sin^ 9] denotes the complete elliptic integral of the first kind, and 
K'{m) = K{1 — m). The complex, PT-invariant potential obtained via an anti-isospectral transformation 
applied to eq. (Q) is 

V^^ix) = -a{a + l)msn\ix + f3,m) - b{b + i)^^!^!i!^±Arn| _ ^2) 

dn [ix + f3, m) 

It is also periodic with a different real period 2K'{m). Furthermore, it is analytically solvable with a finite 
number of band gaps. It is important to understand that the key point for obtaining the above results 
is that unlike trigonometric and other periodic functions, Jacobi elliptic functions are doubly periodic 
functions. This allows both V{x) and V^^{x) to be simultaneously periodic, even though the periods are 
different. Note that the arbitrary nonzero constant /3 in the anti-isospectral transformation, x ^ ix + P, 
is chosen so as to avoid the singularities of Jacobi elliptic functions jTl] . 
Let us first apply our approach to the Lame potentials (b = 0) 

V (x) = a{a + l)msn'^ {x , m) , a = l,2,3, (3) 

which are known to have 2a + 1 eigenstates (band edges) and a band gaps. Let Ej{m) and ipj{x,m) with 
j = 0,1,..., 2a denote the band edge energies and wave functions. The anti-isospectral transformation 
X ^ ix + (3 jlUj yields the PT-invariant potential 

V^'^{x) = -a{a + l)msT?{ix + I3,m) , a = l,2,3, (4) 

with real period 2K'[m). The band-edge eigenvalues and eigenfunctions of V^'^{x) are related to those of 
the Lame potential © by 

Ef'^{m) = -E2a-j{'m) , '4)^'^{x,m) oc tp2a~j{ix + (3,m) , j = 0, 1, 2a . (5) 

Thus, the PT-invariant, periodic potential Q also has precisely a band gaps and 2a -1- 1 band edges at 
energies given by eq. ©. Special mention should be made of the remarkable fact that for any integer a, 
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all bands and band gaps exchange their role as one goes from a Lame potential to its PT-invariant version 

For any band structure problem, an important quantity is the discriminant A jlbj which gives infor- 
mation about the number of band gaps as well as their widths. The question is whether one can relate 
the discriminant for the potential y^^(x) with the discriminant A for the corresponding Lame potential. 
Unfortunately, this is not directly possible by using eq. (0), since it only relates energies of states with 
different numbers of nodes. However, we now derive a remarkable relation using which we can relate the 
two discriminants A and A^^. 

We start from the Schrodinger equation for the Lame potential ^ 

— il)"{x) + a{a + l)msn^(x, m)'ilj{x) = E{m)ilj{x) , (6) 
where a prime denotes a derivative with respect to the argument. On using the relation jl4l I15j 

V^sn(x, m) = —dn[ix + K'[m) + iK{m), 1 — m] , (7) 
and then defining a new variable y = ix + K'{m) + iK{m), the Schrodinger eq. © takes the form 

- ^^J"{y) + a(a + 1)(1 - m)sn'^{y, 1 - m)^p{y) = [a{a + 1) - E{m)]i:{y) , (8) 
so that for the Lame potentials Q we obtain the remarkable relations 

Ej{m) = a{a + 1) — E2a-j{^ — rn) , ^ljj{x^m) oi il}2a-j{ix + K\m) + iK{rn),l — m)^ j = 0, 1, 2a. (9) 
In passing, note that for the special choice m = 1/2, one has several interesting relations: 

Ej{m = 1/2) +E2a-j{m = 1/2) = a{a + I), Ea{m = 1/2) = a{a + I) /2 . (10) 
On combining eqs. © and @ we obtain 

Ef^{m)=Ej{l-m)-a{a + l), j = 0, 1, 2o , (11) 
and hence the corresponding discriminants are related by 

/\P'^{E,m)= l^[E + a{a + l),l-m]. (12) 
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As an illustration, in Figure 1 we plot the real and imaginary parts of the PT-invariant, complex 
potential V^^{x) = — 12msn^(ix + (5,m). Using the well known results for the Lame potential with a = 3 
[T^ and eq. (jSJ, the ground state (lowest band edge) eigenvalue and eigenfunction is 

t/'g(x) = sn(ix + /?, m)[2 + 2m — (53 — 5msn^(ia; + /3, m)] , Eg = —b — — 25^ , (13) 

where ^3 = \/4 — 7m + 'im? . In Table I we have given all the seven band edge eigenvalues and eigenfunc- 
tions. We have subtracted off the ground state energy from the potential so that the lowest band edge by 
construction is at zero energy. Observe from the table that the band edges are both periodic as well as 
anti-periodic with periods 2K'{m) and AK'{m) respectively. 

For the special case of the Lame potential with a = 1, the dispersion relation is also analytically 
known jl7j . We now obtain the dispersion relation for the corresponding PT-invariant potential V^"'^{x) = 
—2msv?{ix + P,m). To that end, we start from the Schrodinger equation: 

- ip"{x) + [1 + m - 2msn'^{ix + (3, m)]ip{x) = ETp{x) , (14) 

where we have subtracted the ground state energy Eg = — 1 — m from the potential so that the new 
potential [F^-^]_(x) = V^'^{x) — Eg has zero ground state energy. On substituting y = ix + [5, eq. (|14j) 
takes the form 

- V'"(y) + [-m + 2msn2(y, m)\il){y) = (1 - S)^(y) . (15) 
Now it is well known that two independent solutions of this equation are given by ^7] 

i?(ix + /3±ai)exp[=F(ix + /3)Z(ai)] 

W{x) = — — — , (16) 

0[ix + p) 

where H, 0, Z are the Jacobi eta, theta and zeta functions, while ai is related to the energy E of eq. (|15|) 

by 

£^ = msn^(ai, m) . (17) 

On using the Bloch condition and the fact that while 9{ix + /3) is a periodic function with period 2K'{m), 
H{ix + 13) is only quasi-periodic ^7j) i-^- 

H{i[x + 2K'{m)] + /3) = H{ix + f3) eiq>[-^K' {m) / K{m)] , (18) 



it is easily shown that the complex, PT- invariant potential [V ]^{x) has a dispersion relation given by 

fc = =F , ±iZ(ai)+i J, (19) 

2K'{m) ^ ' 2K{m) ' ^ ^ 

where ai is given by eq. ((T7|) . 

We now turn to the associated Lame potentials of eq. where without loss of generality we consider 
a > b with both being positive integers. We shall later comment about the case a = b. As has been shown 
by us [121 113j . these are also exactly solvable problems with precisely a band gaps and 2a + 1 band edges. 
However, in many of these cases, some of the bands are unusual in that both the band edges have the same 
period, since some band gaps vanish. In this sense, the associated Lame potentials are much richer than 
the Lame potentials. On using the anti-isospectral transformation, it is easy to see that the band edges of 
the potential ^ and its PT-invariant counterpart 



dn'^{ix + P, m) 



V^^{x) = -a{a + l)msn\ix + /3, m) - b{b + ( , (20) 



are again connected by the relation ((SJ. However, we are unable to relate the discriminants of the two 
potentials since we have not been able to derive an analogue of the relation @. As an illustration, 
let us consider the (a=2, 6=1) associated Lame potential and its corresponding PT-invariant potential 
V^^{x) = —&msv?{ix + P,m) — 2mciP {ix + P , m) / dn^ {ix + P , m) . The ground state eigenvalue and 
eigenfunction is given by 

, , cn(ix + P,m , n, , , ., , , , 

^ X = — ^ 'j;^[3msn^(ix + P,m) -2 + VI^^3rn], E„ = -5 - m - 2^4 - 3m . (21) 

dn{ix + P,m) 

In Table 2 we have given all the band edge eigenvalues and eigenfunctions. 

Finally, let us discuss the associated Lame potentials for the case a = b = integer. In view of the 
well known Landen transformation formula jl4l I18j 



dn(a;, m) + dn{x + K{m),m) = — dn 

a 



X 

— ,m 

a 



1 

a = , , m 

1 + VT^^ 



l-VT 



m 



m 



(22) 



the associated Lame potentials with a = b can be rewritten, apart from an overall constant as Lame 
potentials and so the results derived above for the Lame potentials will go through with a modified modulus 
parameter rh. 



Supersymmetric Partner Potentials: Additional analytically solvable finite band gap potentials can 
be obtained from our previous results by using supersymmetry. The procedure is standard Consider a 
periodic potential V- [x) whose ground state energy is zero, eigenvalues are En "* and eigenfunctions (band 
edges) are ipi, \x)- Let the ground state wave function be denoted by ipg{x) = ip^ \x). One constructs the 
superpotential W{x) = —ip'g{x)/ipg{x). The original potential and its supersymmetric partner potential 
are then given by V±{x) = W'^{x) it W'{x). The eigenvalues are the same for both potentials and their 
un-normalized eigenfunctions are related by 

4"^^(x) ocl/4~^(x) , Vl^)(x)oc ^ + W{x) , n>l. (23) 

This technique is immediately applicable to PT-invariant potentials V^'^{x) of the type given in eq. 1)20(1 
from which the ground state energy has been subtracted. The only caution to keep in mind is that the 
ground state of V^^{x) corresponds to the highest eigenstate of V{x) as indicated by eq. ©. 

Let us apply the above formalism to the Lame potentials. First, for the special case a = 1, one has 

y_(x) = 2msn^(x, m) — m , [V^'^]-{x) = —2'msn'^{ix + P,m) + m + 1 , ■0g(^) = sn(ix + /3, m) , 

W^^ ix) = -I — ^ \ ^ , — [V^ ] + (x) = -2msn^(ix + /? + iK'(m), m) + m + 1. (24) 

sn[ix + fj, m) 

Here, the result of invoking supersymmetry is basically a translation of the independent variable in 
[V^'^]-{x). Such potentials are usually called self-isospectral potentials [T^ . 

For higher a values, the two supersymmetric partner potentials are quite different in shape from each 
other ^5], even though they both have the same band edge eigenvalues. Let us now explicitly consider 
the case a = 3. Here, using the ground state eigenfunction of the PT-invariant potential V^^{x) = 
— 12msn^(ix + I3,m) — Eg, where Eg is given in eq. (|13|) . we find that the corresponding superpotential is 

^ttPT/ \ .cmUx + d,m)dn(ix + B,m) cn( ix + 3, m)sn(ix + 3, m)dn(ix + B,m) 

^ = ~^ — ^ + ^^^^ fTlT^ X R 2r — ' 

sn(zx + p, m) [2 + 2m — 03 — 5msn^(zx + p, m)\ 

so that the supersymmetric partner potential is [y^"^]+(x) = [W^"^(x)]^ + [VF^"^(x)]'. In this way, one 
has discovered another PT-invariant complex potential with a finite number of band gaps. It is plotted in 
Figure 2. 
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We can obtain yet other analytically solvable, complex, PT-invariant potentials by exchanging the or- 
ders of applying anti-isospectral transformations and supersymmetry. For example, we could first determine 
the supersymmetric partner of a solvable associated Lame potential and then compute the corresponding 
PT-invariant potential. The results for the a = 1 Lame potential are: 

= 2msn^(x — K{m),m) — m , [y+]^"^(x) = —2msn'^{ix + P — K{m),m) + m + 1 , (26) 

where the potentials have been adjusted so as to have zero ground state energy. Again, for this a = 1 
example, one finds that is essentially the same as \y^'^\j^{x) with a constant complex shift of 

the independent variable x. 

The situation is much richer and more interesting for Lame potentials with higher a values. For instance, 
it is shown in ref. ^21 that the supersymmetric partner potential of the a = 3 Lame potential is given by 

, , o, , 2m^ sn^fx, m)cn^(x, m) [2m -I- (5i -I- 11 — 15msn^(x, m)l^ 

y+(x) = -12msn2(x,m) + ' 2( \ " R 2/ m2 ' (^7) 

dn (x,m) [Im + Oi + I — bmsii^{x,m)\'^ 

and the corresponding PT-invariant complex potential [V+]^"^ is simply obtained from here by using the 
anti-isospectral transformation x ix + f3 and subtracting off the ground state energy Eg given by eq. H25() 
from it. This potential [1/+]^"'" is plotted in Figure 3. Clearly its band edge energy eigenvalues are simply 
related to those of V^{x) and hence to the a = 3 Lame potential by relation (jSJ. Hence the band edge 
energy eigenvalues of [V+]^"^ are identical to those of and even though the three potentials 

are distinct. For the example under consideration, this is just the statement that the three different complex 
PT-invariant potentials [V^^]^{x), [V^^]+{x), [V+f^{x) plotted in Fi gures 1,2,3 all have the same band 
structure. 

In Table 2 we have given the expression for the band edge eigenvalues and eigenfunctions for the PT- 
invariant complex potential V^^{x) = — 6msn^(ix -|- P,ni) — Imcn^ {ix + P , ni) / dn^ {ix + P , m) suitably 
adjusted by subtracting its ground state energy Eg = —5 — m — 2\/ \ — 3m so that the lowest band edge is 
at zero energy. Using the ground state wave function of this potential, the corresponding superpotential 
turns out to be 

^^,pr sn( ix + P,m)dn(ix + P,m) cn(ix + P,m)sn(ix + P,m) 
W = I z r im- 



cn(zx -|- /3, m) dn(ix + P,m) 

sn(ix + P,rn)dn(ix + P,m)cn(ix + P,m) 

mm — -j^=^= . (28) 

3msn"' (ix + P,m) — 2 + — 3m 



Hence the corresponding supersymmetric partner potential [y^"^]4.(x) = [Ty^"^(x)]^ + [P1^^"^(x)]' is easily 
calculated. On the other hand, yet another PT-invariant potential with the same band edges can be 
obtained by starting from the partner potential of the (a = 2,6 = 1) associated Lame potential V{x) = 
6msn^(x,m) + 2mcn^(x, m)/dn^(a;, m), and applying the anti-isospectral transformation. We get 

[V+f^ix) = -2msn^{ix + /3, m) - 6m ™ ^' "^j + 5 + m + 2V4 - 3m . (29) 

dn [ix + f3, m) 

It is worth noting that again the two potentials and are quite different even though they 

have the same band edge eigenvalues. Further, while the initial associated Lame potential is self-isospectral, 
its PT-transform is not so. In fact, this seems to be true in general. In particular, whereas the associated 
Lame potentials with b = a — 1 are isospectral, we find that the corresponding PT-invariant periodic 
potentials are not self isospectral except when a = 1. 

Finally, let us comment that for the PT-invariant potential i sin^^~''^(x). Bender et al. [H] found that 
the band edge eigenfunctions are always 2tt periodic and unlike other lattice problems, the anti-periodic 
band edge eigenfunctions of period 4tt were absent. They speculated whether this could perhaps be a 
unique signal of PT symmetry. However, in this letter, we have seen many examples where this is not true. 
In particular. Table I shows an example with both periodic as well as anti-periodic band edges, showing 
that the absence of anti-periodic band edges is not a general property of PT-invariant periodic potentials. 

Acknowledgements: One of (US) would like to thank the U.S. Department of Energy and the Interna- 
tional Centre for Theoretical Physics, Trieste, Italy for partial support for this research. 
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Figure Captions 

Figure 1: A plot of the real and imaginary parts of the complex PT- invariant potential [l^^^]_(a;) = 
— 12msn^ {ix+P,m) — Eg, wheie Eg = —5— 5m— 2(^3, and (^3 = V4 — 7m + Am"^ (see ref. ^21)- The potential 
has been defined so as to have zero ground state energy. The plot is for the choice m = 0.75, f3 = 0.5. The 
potential has a period 2K'{0.75) = 3.3715. The continuous curve denotes the real part and the dashed 
curve denotes the imaginary part. 

Figure 2: A plot of the real (continuous curve) and imaginary (dashed curve) parts of the super symmetric 
partner potential [y^^]_|_(x) of the complex PT-invariant potential shown in Figure 1, for the choice 
m = 0.75, (3 = 0.5. 

Figure 3: A plot of the real (continuous curve) and imaginary (dashed curve) parts of the PT-invariant 
potential [y+]^^(x) obtained by first taking the super symmetric partner of the a = 3 Lame potential and 
then applying the anti-isospectral transformation x — > x -|- i(3, for the choice m = 0.75, (3 = 0.5. The 
constant energy —3 — 263, — 26i (see Table 1) has been subtracted off, so that the ground state of [V+]^-^(x) 
has zero energy. 
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Table 1: The eigenvalues and eigenfunctions for the 7 band edges of the PT-invariant Lame potential 
[V^'^]-{x) = -12msn'^{ix+l3,m)-Eg, where Eg = -5-5m-253 [eq. (O], and (5i = y/l-m + Am?, 62 = 
\/4 — m + , S3 = \/4 — 7m + Am'^. The potential has a period 2K'{m). The real periods of various 
eigenfunctions are also tabulated. 



E 


^(-) 




Period 





sn(za; + /3, m) [2 + 2m — (53 - 


5msn^(ix + /3, m)] 


2i^'(m) 


3m + 2^3 - 252 


cn{ix + /3, m) [2 + m — — 


5msn^(ix + (3, m)] 


AK'{m) 


3 + 263- 26i 


dn(ia; + /?, m)[l + 2m — 5i - 


- 5msn^(ix + f3, m)] 


AK'im) 


l + m + 263 


sii{ix + /?, m)cn(ix + /?, m)dn(ix + /3, m) 


2K'{m) 


A63 


sn(ia; + m) [2 + 2m + 83 


— 5msn^(ix + /?)] 


2K'{m) 


3m + 253 + 2S2 


cn{ix + f3,m)[2 + m + 82 — 


5msn^(ix + (3, m)] 


AK'{m) 


3 + 2,53 + 26i 


dn(ia; + /?, m) [1 + 2m + 5i - 


- hm'sr?(ix + /3, m)] 


AK'{m) 



Table 2: The eigenvalues and eigenfunctions for the 5 band edges of the PT-invariant associated Lame 
potential = — 6msn^(ix + /?, m) — 2mcv?{ix + I3,m) /dT?{ix + f3,m) — Eg, where Eg is given 



by eq. (|2H) . Here 64 = \/4 — 5m + m^. The potential has a period 2K'{m). The real periods of various 
eigenfunctions are also tabulated. 

E il^^-^ Period 

dS+ft^) [^"^'"'(^^ + /3, m) - 2 + ^4^3^] 2A''(m 

2^4^:3^- m- 2(54 dngmm) + p^rn) - 2 - m + 6^] 4K'{m 

2V4 - 3m - m + 2(54 dn(S^C) [3^sn^(ia; + /?, m) - 2 - m - (54] 4i^'(m 

4^4^:3^ gS^g [3msn^(i:r + /3, m) - 2 - ^4^3^] 2ir'(m 



5 - 3m + 2^/4 - 3m dn^(ix + /3,m) 2i^'(m 
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